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Abstract 

The effects of meson loops in the vacuum sector of the Nambu-Jona-Lasinio model 
are calculated. Using the effective action formalism we take consistently all next- 
to-leading-order -jy- terms into account. This leads to a symmetry-conserving ap- 
proach, in which all features of spontaneously broken chiral symmetry, such as the 
Goldstone theorem, the Goldberger-Treiman and the Gcll-Mann-Oakcs-Rcnncr re- 
lations are preserved. Contributions to l^q) and F^^ are calculated, and arc shown to 
be substantial, at the level of ~ 30%, consistent with the expansion. The leading 
nonanalytic terms in the chiral expansion of (^g), F^^ and ttIt^ have the same form 
as the one-loop results of chiral perturbation theory. 
PACS numbers: 12.39.-x, 12.39.Fe, 14.40.-n 



1 Introduction 



Over the last few years the Nambu-Jona-Lasinio model [1] has been used 
extensively to study numerous features of strong-interaction physics at low 
energies (for recent reviews on mesons see [2-4], and on baryons see [5,6]). The 
four-fermion point interaction leads to dynamical chiral symmetry breaking 
which plays a key role for the description of low-energy hadron phenomena. 
An attractive feature of the model is that having a small number of free 
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parameters it allows for a unified description of mesons and baryons, with 
both two and three flavors. The resulting hadron spectroscopy has been quite 
successful: the basic mass relations, the electromagnetic form factors and the 
electromagnetic polarizabilities are reproduced reasonably well. 

All the above applications have been performed in the one-quark-loop approx- 
imation. The effective action, the meson and quark self-energies, vertex func- 
tions, etc., are all generated by a quark-loop. This approximation is justified 
by the large- iVc limit, where meson loops are -j^- suppressed. However, meson 

loops cannot always be neglected since in the real world Nc — 3. More impor- 
tantly, the pions are very light, and in many situations pionic effects, although 
formally ;|^-suppressed, are enhanced because of the low pion mass. For ex- 
ample, pionic thermal excitations described by meson loops are dominating in 
hadronic systems at finite temperature. Meson loops are also required for the 
proper description of heavier mesons. For example, the a or the p mesons can 
split into two pions and this effect is again described by pionic loops [7,8]. 

There have been a number of attempts [9-12] to include meson loops in the 
Nambu-Jona-Lasinio model. These approaches are based on expansion in 

terms of Feynman graphs. The symmetry properties of the system, however, 
can be easily violated by an inappropriate choice of diagrams. This is the case 
of Refs. [9-11] where basic symmetry relations, such as e.g. the Goldstone 
theorem, are violated. The reason for this is that the equation determining 
the expectation value of the a field and expressions for the meson propaga- 
tors are not consistent. We show that in order to guarantee the validity of the 
Goldstone theorem on the basis of a consistent expansion it is crucial to de- 
termine both the expectation value of the a field and the meson propagators 
consistently, e.g. using the effective action including one- meson-loop contri- 
butions. Many-body methods which preserve the symmetry properties of the 
theory are called symmetry- conserving approximations [13,14]. The approach 
of Ref. [12] is up to our knowledge the only one which preserves the Goldstone 
theorem at the one meson-loop level. This is achieved by a judicious choice of 
Feynman diagrams. 

In this paper we extend consistently the Nambu-Jona-Lasinio model to next- 
to-leading order in which includes the meson-loop contributions. We use 

the effective action formalism which leads in a natural way to a symmetry con- 
serving approximation. Our approach identifies the Feynman diagrams which 
need to be included to maintain the symmetry properties, at any level of ap- 
proximation. In our calculation the pion remains massless in the chiral limit 
and the basic relations following from Ward-Takahashi identities, such as the 
Gell-Mann-Oaks-Renner relation and the Goldberger-Treiman relation are 
satisfied. Furthermore, the leading nonanalytic terms in the chiral expansion 
of {qq), and ttItt have the correct form, as given by chiral perturbation 
theory at the one- loop level [15]. 

We calculate numerically the contributions of meson loops to {qq) and F^r, and 
find them to be substantial, of the order of 30%. This is consistent with the 
A^c-counting scheme. 

The outline of the paper is as follows: In Sec. 2 we recall the model, estab- 
lish the basic notation and review the formahsm of the effective action. We 
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then proceed to evaluate the effective action at the quark-loop level (Sec. 3) 
and at the meson-loop level (Sec. 4), and derive expressions for the quark 
condensate (Sec. 5), the meson propagators (Sec. 6), and and the pion decay 
constant (Sec. 7) at the one-meson-loop level. In Sec. 8 we verify explicitly 
that our expressions have the correct chiral expansion. Section Sec. 9 presents 
our numerical results. Appendices contain technical details of our calculation. 



2 The effective action 



We use the simplest version of the SU{2) Nambu-Jona-Lasinio model [1] with 
scalar and pseudoscalar couplings only. The Lagrangian of the model is given 

by 



where the point-like four-quark interaction is characterized by the coupling 
constant with the dimension of inverse energy squared, q stands for up 
and down quarks with A^c = 3 colors, and m is the current quark mass. The 
partition function of the system is given by the path integral 



= J Vq^'Vq e-^(«^'«) . (2) 



The integration is over the Grassmann variables q^ and q, and I{q\ q) is the 
Euchdean action 



I" 1 
X{q\ q) = j d^xq^ {dr - iot ■ V + (3m) q - — ({q^(3qf + (qh^-f^rq)'- 



(3) 



The partially bosonized version of the model [16] is obtained by introducing 
auxiliary meson fields $ : 



Z = J Vq^VqV^expl- J d 



1 . 1 ■ 



q^ (a, - ia • V + /3r„$„) q 

(4) 



where we use the following notation 

To = 1; = i75Ta , 

$0 = ^; $a = Pa, a = l,2,3. (5) 
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The quark fields, which appear in the exponent of Eq. (4) in a quadratic form, 
can be integrated out, and the partition function becomes 



Z = e-^ = J D$e-'(*\ (6) 



where / P$ = / VSVP. The so called bosonized Euclidean action I ($) is 
given by 



/($) = -7V,Tr lnD + j d'^x (^^'' - a'm$o + ) . (7) 



where we use the Dirac operator in the form 

D^dr + h (8) 

with the one-particle Dirac Hamiltonian h given by 

/i = -ia • V + f3rA ■ (9) 



The trace in Eq. (7) involves an integration over space-time variables and 
a matrix trace over the spin and flavor degrees of freedom. The trace over 
color gives a factor which we write explicitly. Recalling that the quark 
propagator in the background meson fields $ is given by D~^, we refer to the 
term NcTr In as to the quark-loop contribution. 

For simplicity we use the following notation throughout the paper: the indices 
a,b, . . . contain the field isospin indices a, . . . and the space-time coordinates 
{a,Xa\, b = {P,Xb}, etc. Summation/integration over 
repeated indices is understood. 

Generally, the quark-loop contribution A^^cTr InD can have an imaginary part 
which is related to the anomalous terms. In the case of SU(2) with scalar 
and pseudoscalar mesons this imaginary part vanishes identically and the one- 
quark-loop contribution is given by the real part: A^cTr InD — ^N^Tr ln(£)'^D) 
such that the Euclidean action can also be written in the form 



/($) = -^A^cTr \n{D^D) + J d^x (^^^^ - a^m^o + ym'j . 



(10) 



and DW can be written in the covariant form 

D^D = d^d^ + Z7^r, (9^$,) + $2 . (11) 
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The bosonizcd action /($) is a functional of the fields {x). The effective 
action in the one-meson-loop approximation is given by [17,18]: 

r($) = /($) + -Tr In Vi^- (12) 



The effective action (12) has the following properties [18,14]: 

• At the stationary point of the action, the values of the fields represent 
the ground-state expectation values of the field operators. 

• The inverse field propagators are equal to the second order derivatives 
of the effective action: 

_ g^r ($) 

^'^"-J^^' ^^^^ 

calculated at the stationary point. The same result can be deduced by a resum- 
mation of Feynman graphs, in terms of quark propagators which are dressed 
by a general static potential [14,19]. We refer to the action /($), which is 
leading order in iVc, as the one-quark-loop approximation, and to the effec- 
tive action (12), including both the leading and next to leading order term 

|Tr In as the one- meson-loop approximation. We emphasize that in our 

approach they are both symmetry conserving approximations. 

Most applications of the Nambu ,Jona-Lasinio model so far, have neglected the 
contribution of the second term (the meson loop contribution) to the action 
(12). The purpose of this paper is to include this contribution in the description 
of the vacuum and the pion properties. 



3 The regularized effective action in the one-quark-loop approxi- 
mation 



We first consider the one-quark-loop approximation, in which the second term 
of the effective action (10) is neglected. In this approximation 



r($) = /($) = -^Tr In(L't^) + ^ d^x - a'm$o + ym'^ . 

(14) 

The effective action (14) is ultraviolet divergent and requires regularization. 
To this aim a fermion-loop cut-off Af is introduced. Since the model is non- 
renormalizable, the cut-off Af is kept finite and treated as a parameter. 

In this paper we use two regularization schemes, namely the proper-time reg- 
ularization and the covariant four- momentum 0(4) regularization, as used in 
Ref. [20]. We will compare results obtained with these schemes. The proper- 
time regularized effective action is given by 
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-Tr 



as 



(15) 

In the 0(4) regularization, the four-momentum running in the quark loop 
is hmited by |A;| < Af as explained in App. B.2. In the following we will al- 
ways use the regularized action and skip the index A. The action (15) has a 
translationally invariant stationary point at $ = (Mq, 0, 0, 0) 



57($) 



= 



(16) 



Mo 



The stationary point Mq is identified with the vacuum expectation value of the 
S field in the one quark-loop approximation. Inserting the regularized effective 
action (15) into Eq. (16) we get the following equation for Mq: 



o 



1 - 



m 
Mo 



8iVe5(Mo) =0, 



(17) 



where the function g is defined in App. B. Equation (17) is commonly referred 
to as the gap equation since it determines the energy gap 2Mo between the 
negative- and positive-energy quark states. Obviously, Mo plays the role of 
thge constituent mass of the quark. 

In the one-quark-loop approximation, the inverse meson propagators (13) are 
given by the second variation of the action with respect to the fields, taken at 
the stationary point (Mq, 0, 0, 0): 



Mo 



In momentum space the inverse meson propagators are given by [20] 

K-pi<l') = {4iVe (/(Mo, g2)(g2 + 4MXo) ' 2g{Mo)) + a^} , (19) 



where the function / describing a quark loop with two meson couplings is 
given in App. B. Since we are evaluating the propagator at the stationary 

point (Mq, 0,0,0), we can use the gap equation (17) in order to simplify this 
expression and obtain the pseudoscalar and scalar inverse meson propagators 
in the form 



K-\q')=4NJ{Mo, q')q' + a' — , (20) 

Mq 

771 

K\q')^ANJ{Mo, q'W + ^Ml) + a' — . (21) 
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The propagators have poles at 



— ('21\ 

'"'^ 47VJ(Mo,g2^-m2)Mo ' ^''> 

2 

where m,r and are the on-shell pion and (T-meson masses. In the chiral 
Umit (m 0), the pseudoscalar mesons are massless Goldstone bosons, as 
expected. The physical quark-meson coupling constant 51^^^ is given by the 
residue at the pole of the pion propagator 

glqq- \i^iq' + ml)K^{q% (24) 



which gives in the chiral limit 

y,,, = (47VJ(Mo,0))-t. (25) 



1 



To end this section we recall the iVc-counting rules for quark and meson loops 
in the Nambu-Jona-Lasinio model. From the gap equation (17) it follows that 
the coupling constant ^ must be considered proportional to Each quark- 
loop contributes a factor of A^c, which comes from the trace over color degrees 
of freedom. Thus, the action I ($) is of order A^^- The meson propagators 
contain a quark-loop in the denominator and are hence of order as can be 

seen directly from Eq. (19). From Eq. (25) we find that the physical pion-quark 

_ i_ 

coupling constant gi^gg is of order Nc ^ . The physical pion field (tt = P / g^^qq) is 
obtained by rescaling the meson field P. The propagator for the physical pion 

field has the form KT^{q^) / g^^^, which is of order N^. For each quark-physical 

_i 

pion vertex we have a factor of g^^qq ^ Nc ^ . Analogous relations hold for the 
sigma field. This way, A^c-counting rules in the NJL model agree with those of 
QCD [21,22]. 



The regulcirized effective action including the meson-loop con- 
tribution 



In this section we include the meson-loop contribution to the effective poten- 
tial, which is the second term of the expression (12). This leads to corrections 
which are next-to- leading order in A^c- We now consider the effective action in 
the form (12) 

r ($) = 7 ($) + ^Tr In K'^ ($) , (26) 
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where K ^ ($) is the functional 



_ 5^ 



Evaluating the meson-loop term, |Tr InK ^ ($), we encounter new diver- 
gences which arise due to the integration over the momentum in the meson- 
loop. Let us analyze the meson-loop momentum integrals from a formal point 
of view. In our model, which has only quark dynamical degrees of freedom, 
the meson propagator at the leading iVc-level is given by a chain of quark- 
loops. Wc have regularized the quark-loop contribution introducing a fermionic 
cut-off parameter Af . This regularization, however, does not restrict the four- 
momenta of the mesons. Using the asymptotics of the function / at high 
momenta it can be easily shown that the leading order field propagator has 

the asymptotic expansion K~p{q'^) ^ 6^130^ + (^^^ . Thus, integrating over 

the meson-loop four-momentum in (26) we encounter a quartic divergence. 
Similarly, calculating the pion propagator at the one-meson-loop level we will 
obtain both quadratically and logarithmically divergent terms. We regularize 
the meson-loop integrals introducing a new bosonic cut-off parameter Ab. We 
use a covariant 0(4) regularization for the meson loop, consisting in cutting 
off the meson four- momenta in the loop integrals at — A^. This is the 
simplest possible choice. 

The one-meson-loop effective action (26) has a translationally invariant sta- 
tionary point $ = (M, 0, 0, 0) given by the equation 



The stationary point of the one-meson-loop effective action F is denoted by 
M, in order to distinguish it from the stationary point Mq of the quark-loop 
action defined in Eq. (16). 

We note that M does not correspond to the pole of the quark propagator 
at the one-meson-loop level — it is not the quark mass. It is only the local 
(Hartree) contribution to the quark propagator mass. The quantity M is the 
expectation value of the S field when meson-loop effects included. We may 
rewrite Eq. (28) in the form 



(5r($) 



M 



Sa{M) + -SaUM)K,,{M) = 0, 



(29) 



where we have introduced the following short-hand notation for the n-leg 
quark-loop meson vertices 



/ \ 




Fig. 1. Diagramatic representation of Eq. (31): (a) the quark-loop contribution, (b) 
the meson-loop contribution. The cross denotes the constant term a?{l — m/M). 
The short external lines denote the scalar-isoscalar coupling, the solid line in (a) 
denotes the (regularized) quark propagator, the dashed line in (b) denotes the meson 
propagator K defined in Eq. 33, and the filled blob denotes the meson vertex 5o6c) 
explicitly given in App. A. 

For a — 1,2, 3, the condition (29) is trivially satisfied with a vanishing pion 
field. Prom the variation over the sigma field we obtain 

S'o(M) + ^5o6c(M)i^,e(M) = . (31) 



We represent Eq. (31) digramatically in Fig. 1. The first two terms where 
present at the quark-loop level. The cross is the contact term contribution 
to 5*0 (M), and (a) denotes the quark-loop contribution. The term (b) is the 
meson-loop contribution. The meaning of symbols is explained in the caption. 
Using the iVc-counting rules (see Section 3) it can be verified that the meson- 
loop diagram (b) is of order and hence suppressed by one power in 
compared to the quark-loop diagram, which is of order N^. The leading-order 
meson propagators in Eq. (31) Kbc{M) are evaluated at the stationary point 
of the one-meson-loop effective action, M. In momentum space, similarly to 
Eq. (19), we obtain 



K-^{M, q^) = {4iVe (/(M, q^){q^ + AM^5^o) - 2g{M)) + a^} 

^K-^{M,q') + 5^^A{M), (32) 

where we have defined 



K-^(M,q')^5^^ 



mj{M,q'){q' + AM'S^,) + '''^ 



M 



(33) 



and 



A{M)^-SNMM) + a'{l-^). (34) 
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Since A(M) is exactly equal to the first term So{M) in Eq. (31) we have 

A(M) = So{M) = ~SoUM)KUM) = O (iV°) , (35) 



which shows that A(M) is one order in suppressed compared to the full 
propagators (32) which are of order N^.. Inserting the decomposition (32) into 
Eq. (31) we obtain 



So{M) + ^SobciM)KUM) - ^SoUM)KUM)So{M)K,,{M) + . . . = . 

(36) 

The third term in this expansion is of order and hence one order in 

suppressed compared to the meson-loop contribution (the second term). Fur- 
ther terms are suppressed by even higher powers of Nc- At the present level of 
approximation we should therefore keep only the first two terms in Eq. (36) 
— otherwise consistency in A^c-counting is lost. 

The explicit form of Eq. (36) used in numerical calculations has the form 



2f{M, 0) + -j^ (f(M, Q')(Q' + M- 



X 



K,{M,Q^) 



+ 3 



d 



2/(M,0) + ^/(M,Q2)g2 



0, 



(37) 



5 Quark condensate 



The quark condensate (qq) is given by (qq) = Sr{^)/Sm which from Eq. (26) 
immediately gives 

(qq) ^ -a^{M -m). (38) 



Comparing the result (38) with the corresponding expression at the one-quark- 
loop level, {qq) = — a^(Mo — m), we can see that it is analogous, with Mq 
replaced by M. We shall compare our theoretical results to the empirical 
value [23] 

{uu) = {dd) = -{qq) = -((250 ± 50) MeV)^ . (39) 
2 
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L = 




Fig. 2. Diagrams contributing to the inverse meson propagator in one-meson-loop 
approximation: (a) quark-loop contribution, (b-c) meson-loop contributions. The 
dashed line denotes the meson propagator K, and the filled quark-loop meson ver- 
tices are defined in App. A. 

6 Meson propagators 



The inverse meson propagators, including the one-meson-loop effects, are equal 
to the second order variation (13) of the effective action (26) with respect to 
the fields taken at the stationary point M: 



(40) 



M 



Performing the variation in (40) we can express the propagators in one-meson- 
loop approximation in terms of the functional K (27) and the quark-loop 
meson vertices (30): 



L-,'^K-,\M) - -Sace{M)K,a{M)K,f{M)S,f,{M) 

+\Sabcd{M)K,a{M) . (41) 

The diagrammatic representation of this equation is shown in Fig. 2. The 
first term is the leading-order (quark-loop) contribution K~^[M). We have 
to bare in mind, however, that here it differs from the inverse propagator 
K~^[Mq) in the one-quark-loop approximation, because it is evaluated at the 
new stationary point M of the one-meson-loop effective action. This means 
that the quark propagators appearing in the meson propagator K{M) are 
dressed by solving Eq.(36). The second (b) and third term (c) are the one- 
meson-loop contributions to the inverse meson propagator. The dashed lines 
in the loops correspond to a pion or to a sigma meson. The one-meson-loop 
contributions are suppressed one order in N^, which can be easily shown using 
the A^c-counting rules. 

In Eq. (41) M is the solution of the Eq. (36), in which only the first two terms 
are retained, so as to be consistent with A^c-counting. It is crucial to maintain 
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the same Ai'c-counting consistency in the calculation of the meson propaga- 
tors (41). We include in the inverse propagator L^^ contributions which are of 
order A^^c and N^. In the first term K~^, which is defined in Eq. (32), we keep 
both the contributions of order Nc, and the contribution A(M), of order 
A''^. In the remaining (meson-loop) terms of Eq. (41), the vertices S are of 
order A*";, (because they contain a quark loop) and the leading order of the K 
operators is . Therefore A'c-counting consistency requires the substitution 

of K with K in the two meson-loop terms. Keeping all terms to order we 
obtain the inverse meson propagators in the form 



where the propagators Kah in the meson loops are given by (33). 

Despite the simple form (42) of inverse meson propagators, their numerical 
evaluation for an arbitrary momentum is quite involved. It requires integration 
over the meson- loop four-momentum Q, over the momenta in the quark-loops 
as well as over the proper-time (or Feynman) parameters. Even after tak- 
ing some of the integrals analytically we still end up with many-dimensional 
integrals that must be evaluated numerically. However, we are first of all in- 
terested in pionic properties. Since the pion is very light, we can determine 
the pion mass and the pion-quark coupling constant from the low-momentum 
expansion of the pion propagator. This leads to a radical simplification of the 
numerical calculation. The inverse pion propagator in momentum space can 
be expanded in the form 



Both the constant term, L^^(O), and the g^-term, 2'7r(0), can be expressed as 
a sum of three contributions corresponding to the diagrams in Fig. 2. 

We will now check that the pion remains massless when the meson-loop 

corrections are taken into account. It is crucial for the proof that the expecta- 
tion value of the S field, M, is determined from the stationary-point condition 
for the effective action. We note that at the one-meson-loop level the station- 
ary point of the effective action and the quark self energy do not coincide. In 
the approaches of Refs. [9-11] where the gap equation at the one-meson-loop 
level includes the quark self-energy diagram the Goldstone theorem does not 
hold. It does hold in the work of Ref. [12], where the quark self energy diagram 
has been replaced by a tadpole diagram. Thus contributions of the same order 
in Nc as the calculated meson-loop correction have been excluded in order to 
fulfill the Goldstone theorem. The retained tadpole diagram is exactly the one 
obtained by dressing the quark propagator by a general static potential and 
hence in the approach of Ref. [12] one can recognize a symmetry conserving 
approximation [14]. 




K-,\M) - ^Sace{M)K,a{M)K,f{M)S,f,{M) 

+ ^SabcdiM)K,4M), 



(42) 



(43) 



12 



In order to check the vahdity of the Goldstone theorem we have to evahiate 
the leading-order term L~^(0) in the low-momentum expansion of the inverse 
pion propagator. We write the inverse pion propagator at zero momentum as 
a sum of the contributions of the three diagrams in Fig. 2, in which we have 
now couplings with zero external pion momentum. This leads to a significant 
simplification in the evaluation of the diagrams. The algebra is straightforward 
but tedious, and we refer the reader to Ref. [24] for details. The result is 



N r 

L-\0) = a'-8N,g{M) + ^J d^'Q 



X 



mM,0) + ^,{f{M,Q'){Q' + M')) 



dM^ 



+ 3 



d 



2fiM,0) + ^J{M,Q')Q' 



(44) 



Using Eq. (37) in (44) we obtain immediately 

m 
M' 



(45) 



which is analogous to the one-quark-loop result with Mq replaced with M. 
In the chiral limit (m = 0) we have L~^(0) = 0, showing that the pions are 
massless. Thus, as expected, the Goldstone theorem is satisfied in the one 
meson-loop approximation. This is a general feature. If we derive both the 
gap equation and the meson propagators by performing functional derivatives 
of the effective action, we have a symmetry-conserving approximation. The 
Goldstone pion reflects this. 

For non- vanishing quark current mass m the chiral symmetry is explicitly 
broken and the pion acquires a finite mass. It is given by the position of the 
pole of the pion propagator. Using the inverse pion propagator we obtain 



L-\-ml) = L-\0) - Z^{-ml) ml = 



(46) 



and using Eq. (45) we obtain the on-shell pion mass as a solution of the 
following equation 



2 m 



Since m^r is very small we can approximate Z-,^[—ml) by its value at zero mo- 
mentum. Then we can calculate the pion mass in one meson-loop approxima- 
tion using the results for the low-momentum expansion of the pion propagator 



2 m 
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Q+(l- a)q 



^ ^ / 



Q-aq 



Fig. 3. The routing ambiguity in the meson- loop contributions of Fig. 2, diagram 
(b), to the inverse pion propagator. 

The pion-quark couphng constant can be determined on the one-meson-loop 
level similarly to (24) as the residue of the pion propagator at its pole 



In order to evaluate 771^^ and g^^qq we need, in addition to M which we obtain by 
solving Eq. (36), the wave- function renormalization Z^^iQ). The steps in this 
calculation are as follows: We expand various contributions of Fig. 2 in powers 
of the external momentum q and isolate the coefficient of q^. The momentum 
in loops is cut by the mesonic cut-off A^. In the case of diagram (b) we have, 
for this cut-off prescription, an ambiguity in routing the momentum in the 
loop, as shown in Fig. 3. Depending on the value of the routing parameter 
a, we get different numerical results. In this paper we average with equal 
weights over all a from to 1. This prescription is arbitrary, and one could 
use another one, e.g. choose a = |. As wc will sec, the choice of routing does 
not affect fundamental properties, such as GMOR relation, or the form of the 
leading nonanalytic terms in the chiral expansion (see Sec. 8). Obviously, it 
does affect the numerical value of m^r and to some extent. The algebra and 
final expressions for Zt^{Q) are lengthy, and the reader is referred to Ref. [24] 
for details. 



7 The pion decay constant 



The pion decay constant is an important quantity in our calculations, since 
we use its empirical value in order to fix the parameters of the model. Its 
evaluation allows us also to check explicitly the Goldberger-Treiman relation 
for the quarks, which gives the connection between M, the pion-quark cou- 




(49) 



In the chiral limit we obtain the simple result 




(50) 
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pling constant in the chiral limit. The Gell-Mann-Oakes-Renner relation then 
follows straightforwardly. 

The pion decay constant F.^ is defined by the matrix element for the weak 
pion decay 

miiz) |7r^(g)) = -i^'^^^F.e-^^- , (51) 



where = ?7//75T"/2g is the axial current and tt is the physical pion field tt = 
P/9nqq- Applying the Lehmann-Symanzik-Zimmermann reduction formula we 

can relate the pion decay constant F^^ = F{—rn^) to the form factor F{q'^)^ 
which can be expressed in terms of the expectation value of the time-ordered 
product 

(0|r[.l^(.),7r^(0)]|0) = -iS-^j 70i^(^^)-4^«"^-^- (52) 



To evaluate this matrix element we couple an external axial source in the 
effective action, which gives for the Dirac operator (c/. Eq. (8)) 

L> = /3 {-i^,d, + T^^p + 7^75yJ^) ■ (53) 



Inserting this operator in the action we obtain the generating functional in 
the presence of the source for the external axial field. Then the matrix 

element (52) can be expressed as a variation of the generating functional 



(0|T[A-(z),7r^(0)]|0) = 



1 



5^T 



9nqq J Sj^{z)S7rf^ {x) 



Ln{x) . (54) 



j=0 



The rhs of this equation has the same structure as the one-meson-loop pion 
propagator (40) with one derivative with respect to the meson field replaced 
by a derivative with respect to the axial source. Hence, in the one-meson- 
loop approximation, we obtain contributions to the pion decay constant from 
diagrams analogous to those contributing to the inverse meson propagators 
(Fig. 2), but with one of the external meson couplings replaced by an axial 
coupling. 

We calculate F-^ from the low-momentum expansion as we did for the one 
meson-loop pion propagator, with the same routing prescription. The differ- 
ence is that now we pick up the terms linear in the external momentum g^. 
Also, due to the axial coupling, we obtain different spin-isospin structure of 
the quark-loop functions. The calculation is straightforward but tedious, and 
details are in Ref . [24] . The final result is 

F^ = g^ggMZ^{ml). (55) 
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In the chiral limit this reduces to the Goldberger-Treiman relation [25] for the 
quarks 

g^ggF^^M. (56) 



Using this result and the expressions for the pion mass (48) and the quark 
condensate (38) at the one-meson-loop level we recover also the Gell-Mann- 
Oakes-Renner relation [26]: 

m{qq)^-mlF^ + Oiml). (57) 



We have shown that in our effective action approach the Goldstone theorem, 
the Goldberger-Treiman, and the Gell-Mann-Oakes-Renner relations are valid 
with the meson loops included. Thus the basic relations following from Ward- 
Takahashi identities are satisfied in the one-meson-loop approximation. Again, 
this is a feature of a symmetry- conserving approximation. 



8 Chiral expansions 



In this section we check that the leading nonanalytic terms in chiral expansion 
of the quark condensate, the pion mass and the pion decay constant agree with 
the one-loop results of chiral perturbation theory. 

In chiral perturbation theory [15] we have 



where {qq)o is the value of the quark condensate in the chiral limit. We wish 
to show the same result holds in the NJL model with meson loops. Equation 
(37), which we may write as G{M,m) — 0, defines M as an implicit function 
of m. We have therefore 



(59) 



We find exphcitly 

1^ = a^m/M^ - 8NJg{M)/dM + 0{N^) = ^ + 0{N^) + 0{m) , (60) 

9nqq 



where we have used the identity dg{M)/dM = — 2M/(M, 0), and the relation 
(25). For ^ the only infrared nonanalytic term is the pion- loop term, the last 
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one in Eq. (37). Close to the pion pole this piece is 
-.AT 1 )^.^2^^2 2/(M,0) 



which yields the nonanalytic term ^m'^log{m^/ A^) via Eq. (48). Prom (60) 
and (59) we find immediately 

Using Eq. (56) in (62) we get 



M = M_o l-^^logm^ + ... , (63) 



where Mm=o is the value of M in the chiral limit. From the fact that {qq) = 
—a'^{M — m) we immediately estabhsh the desired result (58). Note, that 
for and F^^ in the above expression we use the leading-A^c pieces. The 
subleading pieces to these quantities are relevant in relative corrections 
to (qq), which is not considered at the present level of approximation. 

Chiral perturbation theory also determines the leading nonanalytic terms in 
the expansion of m,r and F^^. Denoting the leading- iVc values of the pion mass 
and decay constant by m^_o ^-nd F-n-^o (in the notation of Ref. [15] these are M 
and F — the tree-level pion mass and decay constant), we have 

/ 2 \ 

2 2 /i , "^7r,0 



TT 



= F^,o ( 1 - TT-^ <o + - ) ■ (65) 



In order to check that these expansions hold in our treatment of the NJL 
model, we first look at the chiral expansion of ^,^(0). There are cancellations 
occuring between the quark and meson loop contributions, and the final answer 
is 

^+16^^^°^'"''° + -J- ^^^^ 



The above result is independent of the routing ambiguity of Fig. 3. Using ex- 
pansions (63) and (66) in Eqs. (48) and (56) we obtain the desired expressions 
(64-65). 
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We emphasize that the agreement of the results of this section with chiral 
perturbation theory is a direct consequence of the symmetry-conserving ap- 
proximation. 

We have also performed a numerical study in order to check how good is 
expansion (58) in our case. For typical parameters (a = 162MeV, A/ = A;, = 
755McV, ni = 0) we find that (qq) — — (ITS.SMcV)'', and for same parameter 
but m = 13MeV, which fits the pion mass to its physical value, we have 
{qq) = — (175.0MeV)^, i.e. a 5% change. In comparison, Eq. (58) would give 
8% with the chiral scale under the log taken to be IGeV. This shows that the 
leading nonanalytic correction accounts for a major part of the change of the 
condensate as one departs from the chiral limit. 



9 Numerical results and discussion 



At the one-quark-loop level the model has the following parameters: the quark- 
quark coupling constant ^, the quark cut-off Af, and the current quark mass 
m. With the meson loop included we have in addition the meson cut-off Ab. 
We use the empirical values for the observables in the meson sector, namely 
= 93 MeV, and m^r = 139 MeV in order to fix the model parameters. After 
fitting the pion mass and decay constant, we arc left with two free parameters. 
We could fix one of them by fitting the phenomenological value of the quark 
condensate (qq). The results of our calculation for (qq), however, are quite 

sensitive to the cut-off procedure. Furthermore, the experimental range for 

1 

(qq) is very wide, {—^{qq))ixp = 250 ± 50 MeV. For these reasons we do not 
use the value of the quark condensate to fix the parameters. Since we do 
not know the particular physics underlying the regularization of the theory, 
QCD cannot be a guide to fix the loop cut-offs Af and Ab. In this exploratory 
calculation we display results obtained for four values of the ratio Ab/Af = 0, 
0.5, 1, and 1.5. Using the gap equation we eliminate the coupling constant 
1/a^ in favor M, which wc treat as a free parameter. Thus all quantities are 
presented as a function of M for different values of the ratio Ab/Af. 

We present the calculations performed in the chiral limit {m — 0). Introducing 
a nonzero current quark mass m leads to small corrections to the calculated 
quantities. Since at this point wc are mainly interested in the general behavior 
of the quark condensate as a function of the model parameters we find it useful 
to start by considering the chiral limit. The quantity M in one- meson-loop ap- 
proximation is given by the solution of Eq. (31). We use this equation in order 
to determine the fermionic cut-off Af for given values of M and Ab/Af. We plot 
Af as a function of M for different values of the ratio Ab/Af of the bosonic to 
the fermionic cut-off. The results in the chiral limit (m = 0) with proper-time 
and 0(4) fcrmion-loop regularization are presented in Figs. 4 and 5, respec- 
tively. The pion decay constant is fitted to reproduce the experimental 
value. For Ab = we recover the results of the one-quark-loop approximation, 
so that the difference between the lowest curve and the other curves is a mea- 
sure of the meson loop contribution. For Af lower than a critical value we do 

not have spontaneous breaking of the chiral symmetry and Eq. (17) has no 
solution. With increasing mesonic cut-off Ab the critical value of Af increases 
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Fig. 4. The fermionic cut-off Af for proper-time quark-loop regularization as a func- 
tion of M for different values of Ab/Af. The curves correspond to Ab/Af = 0, 0.5, 
1, 1.5 . 
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Fig. 5. The fermionic cut-off Af for 0(4) quark- loop regularization as a function of 
M for different values of Ab/Af. The curves correspond to Ab/Af = 0, 0.5, 1, 1.5 . 
Note that the curve for Ab/Af = 1.5 has the critical Af smaller than for Ab/Af = 1. 

and the curve "moves" up until Ab/Af reaches ~ 1.5 for the proper-time cut- 
off and ~ 1.2 for the 0(4) one. For larger values of Ab/Af the minimal value 
of Af starts decreasing and the curve moves down with increasing Ab/Af. This 
differs from the results of Dmitrasinovic at al. [12] with Pauli-Villars fermion- 
loop regularization where the curves keep moving towards larger values of Af 
as Ab/Af increases (sec Fig. 11 in [12]). For large mesonic cut-offs, Ab/Af over 
2 for the proper-time fermion regularization and Ab/Af over 3 in the 0(4) 
case, we do not have solution of Eq. (36) with fitted to the experimental 
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Fig. 6. The quark condensate y—^{qq)j with proper-time quark-loop regularization 
as a function of M for different values of A^/Af . The curves correspond to A^/Af = 
0, 0.5, 1, 1.5 . The dashed lines mark the empirical bounds. 




Fig. 7. The quark condensate y—^{qq)j ^ with 0(4) quark-loop regularization as a 
function of M for different values of A^/Af. The curves correspond to A^/Af = 0, 
0.5, 1.5, 1 . The dashed lines mark the empirical bounds. 

value. Calculations in the baryon sector of the Nambu-Jona-Lasinio model 
with one quark and zero boson loops [27,28] show that best results for bary- 
onic properties are obtained for M in the range 400-500 MeV, i.e. for values 
of the quark cut-off near the critical value. Although this picture may change 
when meson-loop effects are included, we note that the value M for the critical 
value of the quark cut-off changes very little. For Ab ~ Af the fermion cut-off 
increases by about 30%. 
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Fig. 8. The ratio of the meson- loop to the quark- loop contribution to {qq) as a 
function of M obtained with proper-time fermion-loop regularization. The curves 
correspond to Ab/Af =0.5, 1, 1.5 . 

Comparing the results of the two fermion-loop regularization procedures we 
can see that although the critical values of Af, below which there is no spon- 
taneous chiral symmetry breaking are quite different, the general behavior of 
M as a function of the fermionic and mesonic cut-offs is very similar. 

The quark condensate is given in terms of M and the quark-quark coupling 
constant by Eq. (38). The results in the chiral hmit are shown in Figs. 6 and 7 
for the proper-time and 0(4) fermion-loop regularizations. They are plotted 
as a function of M for the four different ratios A^/Af. The empirical value for 
the quark condensate is obtained in chiral perturbation theory using the value 
of the quark current mass (Gell-Mann-Oakes-Renner relation) [29] or from 
QCD sum-rules. The phenomenological uncertainty (39) is, however, quite 
wide. With both the proper-time and 0(4) cut-offs there is a plateau for M 
between 0.3 and 0.6 GeV and the corresponding value of {qq) depends little on 
Ab/Af. The value in the proper-time case (Fig. 6) is underestimated. This prob- 
lem is known from the one-quark-loop calculations, where better agreement 
can be obtained using a generalized two-parameter proper-time regularization 
function. In the case of the 0(4) quark-loop regularization (Fig. 7) our results 
for the quark condensate are in the phenomenological bounds for values of M 
in the plateau region. This shows that the quark condensate is sensitive to the 
particular quark-loop regularization procedure, a feature well known from the 
calculations in the one-quark-loop approximation. 

In Fig. 8 we plot the ratio of the meson-loop to the quark-loop contribution 
to {qq) as a function of M for the proper-time quark-loop regularization. The 
relative contribution of the meson loop increases with increasing Ab/Af. For 
values of the parameters in the range 0.3 < M < 0.6 GeV and Ab/Af < 1 we 
have (^?)meson/(??)quark < 40%, 1.6. of the ordcr of as expected. 

Looking at Figs. 6 and 8 we see that as Ab/Af increases the meson- loop con- 



21 



0.8 
0.7 

^ 0.6 
I 0.5 
^0.4 
, I . 3 
0.2 
0.1 


0.2 0.4 0.6 0.8 1 

M [GeV] 

Fig. 9. The ratio of the meson-loop to the quark-loop contribution to as a function 
of M obtained with proper-time fermion-loop regularization. The curves correspond 
to Ab/Af =0.5, 1, 1.5 . 

tribution to the quark condensate increases, while the quark-loop contribution 

decreases. The sum of the contributions remains roughly constant. Eq. (38) 
shows that the inclusion of the meson-loops does not significantly change the 
four-quark coupling constant ^. 

Let us now come back to the A^c-counting scheme, which we extensively used 
throughout the meson-loop calculations presented above. We have calculated 

the meson-loop corrections to the quark condensate, the pion mass and the 
pion decay constant, showing that they are one order in suppressed. How- 
ever, the expansion in powers of should be used with care since N^. = 3. 

Furthermore, the meson-loop corrections depend on the mesonic cut-off 
and are infinite for Ab — > oo. Hence it is important to check how well does the 
iVc-counting scheme work for the particular values of the parameters. 

In the previous section we have compared the meson-loop and quark-loop con- 
tributions to (qq) and found that for reasonable values of the parameters the 
ratio (gg) meson/ (^Q') quark IS in the expected range. Here we compare the meson- 
and quark-loop contributions to F,^. Since we have fixed the parameters in or- 
der to reproduce the experimental value of the pion decay constant the sum of 
both contributions gives F7r,quark + -^7r,meson = 93 MeV, but their relative con- 
tributions depend on the choice of parameters. The meson-loop contribution 
is negative. In Fig. 9 we show the ratio |F7r,meson/-p7r,quark| with proper-time 
quark-loop regularization. The relative contribution of the meson-loop grows 
with increasing M. It is larger for larger values of Ab/Af. As we mentioned 
in the previous section, the empirical value of F^^ can be only reproduced for 
Ab/Af less than about 2, which sets an upper bound for the meson-loop cut-off. 
For values of Ab in this range and M < 0.6 GeV the meson-loop correction 
docs not exceed 40% of the A^^c-leading quark-loop term. This is important 
because it keeps the corrections in the expected range. It justifies a pos- 
teriori the use of the A^c-counting scheme. The above observations are also 
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true for the quark condensate (see Fig. 8) as well as for the 0(4) fermionic 
regularization. In the case when Ab — Af and 0.3 < M < 0.6 MeV we find 
that the meson-loop contributions are of the order oi ^ = ^ compared to the 
leading-order contributions. 



10 Conclusion 



We have shown that the effective action formalism in the Nambu-Jona-Lasinio 
model leads to a symmetry-conserving approximation, allowing us include 

consistently meson-loop effects while preserving the usual properties associ- 
ated with spontaneous chiral symmetry breaking (the Goldstone theorem, the 
Gell-Mann-Oaks-Renner relation, the form of the leading nonanalytic terms 
in the chiral expansion of various quantities). Indeed, meson-loop effects will 
destroy these properties unless both the gap equation and meson propaga- 
tors arc treated consistently. This has already been noticed in Ref. [12], where 
Feynman diagrams were used. In our approach the conservation of the sym- 
metry properties is a natural consequence of the consistent application of the 
expansion. We have found that though meson-loop contributions are 

suppressed as compared to the leading-order quark-loop contributions they 
lead to substantial corrections. For the physically reasonable values of the 
model parameters the effects of meson loops to the pion decay constant and 
to the quark condensate are of the expected order of 30%. 
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A Quark-loop meson vertices 



In our meson-loop calculations we need the meson vertices which couple dif- 
ferent number of mesons (we need up to four) through a quark loop. They 
are given by the variations of the effective action with respect to the meson 
fields. To obtain these we first expand the action in powers of the meson field 
fluctuations around the stationary point. We perform the expansion for the 
proper-time regularized effective action. Using the proper-time expressions we 
can easily obtain the four-momentum regularized expressions as well. 
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We start from the proper-time regularized effective action 

AT °° J 2 



wfiere D^D is given by 

D^D = d^d^ + Z7^r« J + $2 . (A.2) 



The fields $ are now fluctuating around the stationary-point constant fields 
which are of the form ^^"^ = {S, 0, 0, 0} 

$ = + 5$ . (A.3) 



Expanding D^D in terms of the meson field fluctuations we obtain 

D^D = G^' + V, (A.4) 

where 

G-i = a^a^ + (A.5) 



is the zero-order term in (5$ which is diagonal in momentum space and the 
fluctuation term is given by 



y« = 2$:*(5<i>, + i7^r«(a^5$,) , 

y(') = (5$)' (A.6) 

and consists of terms of first, V^^\ and second order, V^'^\ in the field fluc- 
tuations. Next we expand the exponent in the fermionic part of the effective 
action (A.l) in powers of V using the Feynman-Schwinger-Dyson formula and 
the cyclic property of the trace 
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.3 1 



-jjdu J (i^;l/e-*(^-"-'')^o Ve-^"^oVe-'^^o' (A.7) 



^ 1 1— M l — U — V 

+^J(iuJ dv J (i^/;ye-'(^-"-''-"')^o Ve-'"^o Ve-'^^o Ve-*"'^o' 



The first term in this expansion (zero order in V) gives the quark contribution 

to the one- quark-loop effective action (15). The term hnear in 5$ cancels with 
the mesonic terms in the expansion of the action since we expand around the 
stationary point. 

The term quadratic in 5$ acquires contributions from the linear and quadratic 
in V terms as well as from the mesonic part of the effective action. Evaluating 
it in momentum space and taking the second variation with respect to the 
meson fields we obtain 



S - 



(5$„(gi)5$^(g2) 

= SapS{q^ + q^) {47Ve [f{S, q,){ql + 4S'6^o) - MS)] + a'} , (A.8) 
where the functions / and g are calculated in App. B. 

The third order term in the expansion of the action / in powers of the field 
fluctuations 6^ acquires contributions from the second and third order in V 
terms in (A.7) 



2 CO 1 

/(3^) = iV,^Tr J cis||dKl^We-^(i-")^oV(2)e— «o\ (A.IO) 



Cv 



A-2 

oo 2 ^ 1— « 



A-2 



(A.11) 



where ^ denotes the sum over the permutations of F^^^ and V^'^\ Varying 

Cv 

this expressions with respect to the meson flelds we obtain the corresponding 
three-leg quark-meson vertex functions in momentum space 



(A\ 53/(3^) 



25 



A 



X J dss|d«|-|^e-^(i-'^)«'=+^i)'+^')e-^^(^'+^'\ (A.12) 



X (4$^„*$J$^* + 5ap^'^ qi-q2 + Sp^^'^ g2 • ?3 + (^a7^? ?1 " 



X dss^ j du j dv j 



(27r)4 

A-2 

^g-s(l-u-^)((fc+qi)2+52)g-s«((fc+gi+g2)2+52)g-s^(fc2+52)^ ^^-^3^ 

where we have introduced the notation ^ for the sum over the permutations 

p 

of the couples of meson-field indices and momenta {{a, qi), (/?, ^2), (7, qa), ■ ■ ■}■ 
Summing the two contributions we obtain finally 

Sai>c^sitl+Sif^ . (A.14) 



The fourth order term in the expansion of / in powers of the field fluctuations 
5$ acquires contributions from the second, third and fourth order in V terms 
in (A.7) 



/(4)^j(4.l)+j(4S)^^(4C)^ (A. 15) 

00 1 

I^^^^^N.Tr I ds|ydiiy(2)e-^(^-")^oV(2)e-^"<, (A.16) 



A- 



-.2 1 



li^B) = _iv, ^ Tr J dsjl du j dv 

Cv ,-2 

^y{'2)^-s{l-u-v)G-^y{l)^-suG-^y{l)^-svGQ^ ^ (A.17) 



00 3 -'^ 1— n 1—u—v 



j(4C) ^ ^j.^ J dsj J du J dv J 



dw 



A-2 



x"\/(i)e"*(^~"~''~'^)'^o V(^^e"^"'^o V^^^e"^'"^o V^^^e"'''"'^o \ (A. 18) 

We evaluate the traces in momentum space and vary with respect to the 
meson fields in order to obtain the corresponding contributions to the four-leg 
quark-meson vertex 
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^4j(4A) 4! 

'S'iid = / \r>F. / / \ / \ = V 5(^1 + ^2 + ^3 + ^4) 

(5$«(?i)5$^(g2)(5^7(?3)(^^5(?4) p 



x5,^5^5| dss J du J -^^e-^^'--'>'^^'+''^+''^^^^^^^ (A.19) 



(B) 54j(4B) 4! 



xSafiSysiqs ■ q4: + 4:S'^5^o) j dss^ j du j dv J 



00 ^ ^ 

"f 

^^-s{l-u-v){{k+qi+q2)'^+S^)^-su{{k+qi+q2+q3)'^+S^)^-sv{k'^+S^) ^ (A.20) 



^4^(4C) 4! 

r ^ 

Cyi 

x((?i • 52) (?3 • ^4) - {qi ■ g3)(g2 • 54) + (gi • 54) (g2 • 53)) 

00 1 1— « 1— «— 1' ^ 
X J dss'Jdu J dv J ci^^|Ale-^(^-"-^--)«'=+«^)'+^') 



(27r)^ 

^g-s«((A;+gi+g2)^+S^)g-st;((fc+gi+g2+93)^+S'2)g-s-i<;(fc2+S'2) ^ (A.21) 



A-2 



6 

where in the last contribution ^ stays for the sum over combinations, origi- 

Cyl 

nating from different orderings of the two terms in V^^^ (A. 6) 



X:(<^^.$J$^* qi ■ qj) = Sap^^^f qi-q2 + S^,^$^f qi ■ q^ (A.22) 

Cyl 

gi • 54 + 5/37 *a*? ?2 • 53 + 5/35$a $7 ^2 " ^4 + 5^5*a ^3 " g4. 

Summing up all contributions we obtain finally 

Sabcd — 'S'abcd + 'S'lftcd + ^ibcd ■ (A. 23) 
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B RegulEirization functions 



Here we calculate the regularization functions for the quark-loop vertices in the 
diagrams for the gap equation and pion propagator. We start by evaluating 
the regulators in the proper-time scheme. The 0(4) regulators can be then 
easily obtained using the intermediate results of the proper-time case. 



B.l Proper-time regularization 



We start by evaluating the function g emerging in the one-quark-loop gap 
equation. It is given by 



' ^A^,e-^-S^EJ^]], (B.l) 



167r2 V V^f 
where Ei is the exponential integral defined by 



°° -xt 



E^(x) ^ J dt^-^ . (B.2) 
1 



The function / is defined by 



(27r) 



A-2 



^ * 



S' + «(l-«k!V (B.3) 



B.2 0(4) regularization 



Here we define the quark-loop 0(4) regularization. The running four-momen- 
tum of the quark loop is limited in the following manner. We obtain the quark- 
loop vertex functions in our O (4) regularization scheme using the intermediate 
results from the proper-time calculation. After having completed the square 
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for the fermion-loop four-momentum k in the proper-time expressions we take 
the hmit in the proper-time fermionic cut-off A^"^ oo using the identity 

oo 

/n' 
ds s"e-"^ = — -— , forn > 0, ^ > . (B.4) 
zln+l ' — ' ^ ' 



(AD" 



Then we cut off the integral at some 0(4) cut-off Af^^^^ and obtain the 0(4) 
regularization functions. Note that in the following expressions both the k'^- 
and M-integrals can be taken analytically. This leads, however, to quite lengthy 
sums of logarithms and rational functions, and for the sake of clarity we present 
the regularization functions keeping the integrals. 

The function g{S) can be easily obtained from the 0(4)-regularized action 

A2 




Using Eq. (B.3) we obtain 

A? 1 
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